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' We analyze the itinerant model for antiferromagnetism, which was developed previously by Plis- 

chke, Mattis, Brouers and Mizia. In this model we include both; single-site and two-site electron 
■ correlations. Including additionally band degeneration into the model allows for considering intra- 

' atomic exchange interactions in the Hamiltonian. The modified Hartree-Fock approximation for the 

two-site interactions is used. This approximation gives the spin band narrowing, which is the same 
for both spin directions and affects possibility of antiferromagnetic ordering. We use the Green 
function technique and CPA decoupling. This allows for the change in shape of the spin bands, 
described by the correlation factors, which decreases kinetic energy of the system. The effective 
Hartree field and the gain in kinetic energy due to the on-site and inter-site correlation factors drive 
, the antiferromagnetism. 

s. 

^ 1. Introduction 

a , 

The basic model for magnetic ordering of itinerant electrons in solids is the Hubbard model [1]. In the mean- field 
approximation the Hubbard model leads to the well-known Stoner model for magnetism. Extension of the Hubbard 
model is the model, which in addition to the Coulomb repulsion C/ includes also Coulomb-" exchange" interaction 
i J = | 1/r | jy,iX). The solutions of the 3-dimensional Hubbard model did not yield the exchange field small 

enough to be justified physically, therefore to obtain the magnetic state Hirsch ( [2,3]) and others (see e.g. [4]) included 
also the nearest-neighbor inter-site interactions; (i, j \ 1/r | k, I), where (k, I) = and i, j are the nearest-neighbor 
lattice sites. In our model we introduce into Hamiltonian the on-site Coulomb repulsion Uo single-site two sub-band 
l/^ i interactions and also the inter-site nearest-neighbor interactions. Analyzing the inter-site correlations we deal with 
probabilities defined for the product of four operators, in which two-operator product is replaced by its stochastic 
value. We have interpreted the stochastic value of two-operator product as the product of probabilities for electron 
transfer between two neighboring atoms, which is the product of probability that site i is empty for electron witch spin 
+cr; 1 — n a , and the probability of site j being occupied; n a . This interpretation makes the connection of standard 
probabilities used in the CPA method with Hirsch's average bond energy for spin <r; I a = n a {l — nf). The connection 
allows us to expand our model to more realistic cases of the weak inter-site interaction being accompanied by the 
weak or strong single-site Coulomb interaction. 
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2. Hamiltonian and the Coherent Potential Approximation 
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The Hamiltonian for one degenerate band can be written in the form, which we already previously introduced (see 



[5]) 



H = - ^2 Uj (cf X(r c jva + h.c) - fi ^2 *W + ( iX >3 l/ \ 1 / r \ kl J> lri7 ) c i\a c tvv' Cl w' Ck ti°> W 



a, a' 



where tj» is the nearest neighbors hopping integral, /io is the chemical potential, cf XlT (ci\ a ) creates (destroys) an 
electron of spin a in a Wannier orbital A on the i-th lattice site, the indices A, v, w, 1? numerate the sub-bands in the 
degenerated single band. Taking into account in Hamiltonian (1) only single-site and two-site interactions, and also 
single sub-band and two sub-band interactions we obtain and retain the following matrix elements; 
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• single-site, single sub-band interaction 
Uq = (iX, i\\l/r\i\, iX) , 

• single-site (subscript "in"), two sub-band interactions (for A ^ v ) 

V in = {i\,iv\\/r\iX,iv), J in = (iX,iv\l/r\iv,iX), J' in = {iX,iX\\/r\iv,iv) , 

• two-site interactions 

Vo = (iX,jv\l/r\iX,jv), J = (iX,jv\l/r\jv,iX), 
Jq = (iA,iA|l/r|jV,jV), At = {iX,iX\l/r\ju,iX) . 

Keeping these interactions in Hamiltonian (1) we obtain 



H = -t ^2 {cjgCja + h.c.) - no rij + U y~] n ia ni- a + V mnj 

<ij>a i ia <ij> 

+ ^ At {4a c j<r + h - C ) ( n i-v + n i-a) + J C ia C ta> C ™' C 3<T + J ' YL C ia C t-cr C 3-a C 30i ( 2 ) 

<ij>a <ij>aa' <ij>a 

where in the degenerated band model in the Hartree-Fock approximation we have 

U = U a + ( P -l)(J rn + J / m + V m ), (3) 

and p is the number of degenerated sub-bands in the band. The effective inter-atomic constants J, J', V, At for 
the degenerate band are given by pJo, pJ' , pVo, pAt respectively. We will apply mean field approximation to the 
inter-atomic weak interactions J, J', V and At. After these approximations the Hamiltonian (2) will take on the 
following form 

h = - % ( c t c 3° + h - c ) - Y2 ^° n * + Y2 M<Tn ^ + U Y1 ™^ n *-^ ( 4 ) 

ij i ia ia 

a 

where the mean (molecular) field is given by 

M a = JhJ* + Vn + At (c+ a Cj- a ) (5) 

and the spin dependent hopping integral tfj is expressed as 

tl = t - At (n,_ CT + ) + ( J + J') I~r° + ( J - V) 1% , (6) 

where h i(7 is the average number of electrons with spin a on thei-th lattice site, and the parameter If- is defined as 
^ij = ( c tcr c j<?)- 

In the case of ferromagnetism parameter If- was spin dependent (see e.g. [2,5]). For the antiferromagnetism indices 
i, j belong to the neighboring sublattices a, (3 with opposite magnetic moments. In the result If- = If- a = Iaf is spin 
independent. 



3. Application to antiferromagnetism 

For the antiferromagnetism we will use the diagonalization of Plischke and Mattis [6], Brouers [7] and Mizia [8]. 
The crystal lattice will be divided into two interpenetrating sublattices a, (3, with the average electron numbers equal 
to 

n- = ^, »r = ^, (7) 

where magnetization m — n a a — n~ CT = nj^ — n^. We use mean field approximation for the interactions between the 
neighboring lattice sites. For the stronger on-site interaction we use the coherent potential approximation (CPA). To 
describe the antiferromagnetic state we write the model Hamiltonian (2) taking into account the two interpenetrating 
magnetic sublattices , a, (i, with magnetic moments given by Eq. (7) 
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H = - E E *i3 ( c icr c J° + h - c ) ~ E (*" ^ M a ) n irT + ^2lln ia n^ a + ^Un la n^ a . 



(8) 



Next we use the main idea of the CPA formalism [9] and split the above stochastic Hamiltonian into the homogeneous 
part 



Ho = - E E 1 % ( c *W + h - c ) -E^^ + E^^ + E E ^ n - 

iEa jef3 ic i£a i£fj 



and the stochastic part 



ie/3 



(9) 



(10) 



where t^ CT is the stochastic potential given by 

V a = 



(11) 



M = Mo — ^ n — an d S£ are the self-energies on sites a; = a(/3) for electrons with spin a. Using Hamiltonian 

(9) in the equations of motion for Green functions we obtain the following equations in the momentum space 



s + [i - Sk 
-s k e + fj, - 



G^(k,e) G^(k,e) 



where 1 is the unity matrix. 

Next we use these Green functions in the relations for electron numbers 



oo 



(12) 



where for x = a (/?) we have respectively 



' 11 [ 11 ^\G (k,s eff )-G°(k,-e eff ) 



2 y e + fj, - £ 



If we use these results and changing summation over "k" to integration, and inserting for self-energies their first order 
approximation 



X aW =F AF ^ T F AF ^=F AF ^ T A, (13) 



we obtain 



oo 

i±° = \ J p {e)^[S ± "f{E-n) + S^f{-E-n)]dE, 



where E = \fe 2 + A 2 , S ±a = y and /(-E) is the Fermi distribution function. The antiferromagnetic energy 

gap A is given by the expression; A = F AF ^, in which F AF is the value of the effective exchange interaction 
F AF = U — zJ. It can be seen that contrary to the case of ferromagnetism (see [5]) the inter-atomic exchange 
interaction, J, in effective field F AF is opposing the antiferromagnetism. The antiferromagnetic moment per atom 
(in Bohr's magnetons) will be given now by the following expression 
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m = ^ E p r) lf(Ek) f(~E k )], (14) 

k 

where P^ a is the occupation probability of state {k, ±a), Ek = ^Je'f. + A 2 , and 

P£° = (15) 

The effective electron dispersion relation e^ a used above has the following form 

e^^s k S^, (16) 

where 



Jl J , C±<T _ . i E k T A 



e fc = e k L A F (n,m) , 5 fc CT = 



S fe ± A' 



Laf (n, m) = 1 - n— + — (-2 J — J' + V") (m) . 
to to 



The chemical potential \i may be determined from a carrier concentration n on the base of the equation 

k 

Inter-atomic interactions J, J', V and At used in this model change the effective half bandwidth according to the 
relation, where D = D LAF( n , m) is a factor of the bandwidth change. For the antiferromagnetic ordering the change 
of width of the band is independent from spin direction, since 

Iaf (m) = 1 J^~ k if(Ek) f(-E k )\. (18) 

k 

Inserting Eq. (16) and Eq. (15) to Eq. (14) we obtain 

1 = - pAF ^ E [f(Ek) f{ - Ek)] - (19) 

k 



4. Numerical examples 

We calculated the critical values of interaction F£F from Eq. (19) for different values of assisted hopping At 
assuming that the other inter-site interactions are equal zero. We used constant density of states. The assisted 
hopping a = At/t corresponds to the solid line in Fig. 1. Different curves calculated from complete Eq. (19) for 
different values of assisted hopping a are also shown in Fig. 1. One can see that with growing a the critical exchange 
constant necessary for AF is decreasing for small hole concentration. The decrease is even stronger in the split band 
limit at n — > 1 . This effect could possibly explain the existence of the AF at the almost full subband in the high 
temperature superconductors. There is one limitation in here; that these are s-wave (so + s x 2 +y 2^ superconductors, 
since the assisted hopping contributes only to the superconductivity of this type of symmetry (see [10] for more 
detailed analysis of this problem) . 



4 



0.25 0.5 0.75 1 1.25 1.5 1.75 2 
n 



FIG. 1 



FIG. 1. Dependence of critical field on carrier concentration F cr for constant density of states with assisted hopping only; 
a — - solid line, a = 0.3 - dashed line, a = 0.4 - dotted line. 
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